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Abstract 

We present and numerically solve a model of the photorefractive Gunn effect. 
We find that high field domains can be triggered by phase-locked interference 
fringes, as it has been recently predicted on the basis of linear stability con- 
siderations. Since the Gunn effect is intrinsically nonlinear, we find that such 
considerations give at best order-of-magnitude estimations of the parameters 
critical to the photorefractive Gunn effect. The response of the system is much 
more complex including multiple wave shedding from the injecting contact, 
wave suppression and chaos with spatial structure. 
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I. INTRODUCTION 



Segev, Collings and Abraham (SCA) |1J proposed an interesting mechanism for producing 
Gunn domains by means of a photorefractive parametric excitation. It occurs when two 
optical waves of slightly different frequencies are incident upon a biased semiconductor 
crystal doped with deep impurity centers. The authors conjecture that the resulting traveling 
interference pattern excites multiple high-field Gunn domains which move phase locked 
with the interference fringes. Recently, L. Subacius, et al, |2|] proposed an efficient way of 
creating simultaneously a number of quasilocalized high field Gunn domains through hot 
carrier transport in spatially modulated and nonuniformly heated electron-hole plasma, and 
present some numerical results and preliminary experimental confirmation. 

In this paper we present a consistent model of the photorefractive Gunn effect and carry 
out numerical simulations to understand the dynamics of the system. We find that high 
field domains can indeed be triggered by phase-locked interference fringes, as suggested by 
SCA. However, the response of the system can be very complex, and it is not possible to use 
a simplified version of Kroemer's NL criterion, as suggested by SCA, to predict the number 
of high field Gunn domains traveling through the sample. Indeed, our results indicate that 
with appropriate values of the parameters of the system the response becomes chaotic and 
this is, therefore, another example of driven (sinusoidal interference pattern of intensity I) 
chaos. 



II. MODEL EQUATIONS 

The following equations describe the photorefractive Gunn effect, 

SI(N D - Ni) - 7 nN l D , (1) 

(2) 



dN* D 

dn _ dNh \dJ_ 

dt dt q dz : 



f = -± (n + N A -K) (3) 

where (|l|-|2[) are the continuity equations for ionized donors and for electrons respectively, 
and (§) is Poisson's law. In these equations, z, is the space variable in the direction of current 
flow, t, is the time, N l D (z,t), represents the number density of ionized donors, N D , the total 
donor density, I(z,t), the incident light intensity, S, the photoionization cross section, 7, 
the recombination rate, n(z,t), the electron number density, — q the electron charge, J(z,t), 
the current density, E(z,t), the space-charge field inside the crystal, Na, the density of 
negatively charged acceptors, and, e s , the low frequency dielectric constant. 
The light intensity, I(z,t), is given by 

I(z,t) = I [1 + m cos(Kz + Qt)} (4) 

which, as decribed in 0, is the intensity of the moving interference pattern formed when two 
quasimonochromatic plane waves of slightly different frequencies, uo and u + Q (Q << ui) 
and slightly different angles of incidence illuminate a bulk semiconductor crystal. In (f|) 
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K = 2tt/A is the interference wave number, m the modulation depth of the interference 
grating, and Jo the total average intensity. 

Following Sze || [Eq. (28) in chapter 11] the current density, J, includes drift and 
diffusion terms, 

d[D(E)n] 

J = qnv(E) +q 1 j- J J . (5) 

which is the standard form of the drift- diffusion current density [^H where, v(E), is the 
electron drift velocity, and D(E) the diffusion coefficient. The drift velocity of the electrons 
is a known function of the electric field exhibiting negative differential resistance. In the 
following analysis we use a saturating drift velocity function given by, 



v(E) 



E/E s 



A(E/E S 



(6) 



where v s is the saturation drift velocity, E s the saturation field, and A and f3 dimensionless 
constants that depend explicitly on the mobility of the material (fi = v s /E s ). 

Differentiating the Poisson equation @ with respect to time, inserting the result in 
the continuity equation for electrons (H), and integrating with respect to space, results in 
e s dE/dt + J = Jtotai- Here the constant of integration, J total-, is the total current density. 
Introducing the electron current density given by (|5]) in the previous equation, leads to 
Ampere's equation, 

d[D{E)n] 8E 
Jtotai = qnv(E) +q g + e a ^. (7) 

Notice that in Eq. (6) of |TJ the current density erroneously includes the displacement term, 
and this leads to an incorrect Ampere's equation in which the displacement current drops 
out. 

In addition to these equations the electric field distribution must satisfy the reverse bias 
condition for a given applied voltage V, 

L 

Edz = V. (8) 

o 

Thus, to model the photorefractive Gunn effect we use the correct form of Ampere's law 
(^), the continuity equation for ionized donors (|TJ) , Poisson's law (||), the bias condition 
(§), and appropriate initial and boundary conditions. Solution of these equations using the 
known functions for the light intensity (||) and drift velocity provide the four unknowns, 
N l D (z,t), n(z,t), E(z,t) and J to tai(t). 

It should be emphasized that boundary conditions play an essential role in the existence 
of the Gunn effect. Although it can be debated which are the correct conditions to apply 
in order to simulate a particular experiment, periodic boundary conditions, as used in [|TJ, 
should be avoided: they yield a total current density which is constant in time when a dc 
voltage bias is imposed. In fact, integrating Ampere's law we find Jtotai = L~ x \e s dV / dt + 
Jq Jdx) . Then imposing periodic boundary conditions and a dc voltage bias we find that 
d Jtotai /dt = (c/ L) ^{dJ/dx)dx = (provided we have a domain moving at constant speed 
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c; note that with physically reasonable boundary conditions, the domain cannot move at 
constant speed as it arrives at a realistic contact, e.g. Ohmic). Since the Gunn effect refers to 
time-dependent oscillations of the current under dc voltage bias, one should not use periodic 
boundary conditions when discussing it. 

The most important boundary condition in a study of the Gunn effect in long samples is 
that for the injecting contact. In fact, the formation of Gunn domains is due to a periodic 
destabilization of the boundary layer attached to such contact during the oscillations . 
During the formation of a new wave at the injecting boundary the displacement current 
plays a crucial role and can not be neglected. Thus, we use an ohmic condition at the 
injecting contact, E = pJ to tai [@HH' anc ^ dE/dz = at the receiving contact (which is thus 
passive and integration time is saved). When p is such that E/(pqN A ) intersects v(E) on 
its decreasing branch, the Gunn effect is found for m = and appropriate values of the bias 

v IB. 

In order to solve numerically our model equations, we will first write them in nondimen- 
sional form. Let us redefine our variables in dimensionless form as follows: 



1 v s „ E 



(9) 



J total N l D n v{FE s ) 

J = tt, P = tt> = — , v{F) = . 

v s qN A N A N A v s 

Inserting this in (1), (3)-(8), we obtain 

£ (f + '!)=•««>-"» *<n ("» 

dF 
dp 

— = io [1 + m cos (ky + ws)] (1 — ap) — fir/p, (12) 
i 

^ Fdy = 0. (13) 

Here we have defined the following nondimensional parameters (we assume constant diffu- 
sivity, D(F) = D): 

e s E s ± 1 qDN A e s E s 

V, 5 = FT; P = 7, 



qN A V E S L ' v s e s E s ' ev s q 

N A . Se s E s N D e s E s e s E s 

a = , to = ^--'05 w — k = 

N D ev 8 qN% ev s qN A eqN A 



(14) 



III. RESULTS 

We have solved (|T0D - (|1^) numerically for different values of the parameters w and 
i , keeping fixed e = 4.138 10" 4 , m = 0.1, p = 1.8, k = 2vr, 5 = 0.05, a = 0.01, (3 = 
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3.5294, and = 3. These numerical values correspond to those used by SCA jj]]. We were 
particularly interested in testing a central piece of SCA's analysis, namely their particular 
use of Kroemer's NL criterion, Eq. (20). Prompted by SCA's suggestion that the length in 
Kroemer's criterion may be the distance between multiple domains, we combine Eq. (20) 
with Eq. (17) for the electron density n\, and use I = L/N, where £ is the distance between 
adjacent domains. Then we obtain a formula for the maximum number (N) of high field 
domains which may coexist: 



We have used SCA's numerical values (Jo and Q are measured in W.m -2 and Hz, respec- 
tively). Fig. 1 depicts the level curves Q(Iq) which are obtained when N Q takes on different 
integer values. If SCA's theory holds, N = i between the level curves N = i — l and N = i. 
However, we observed Gunn domains (Fig. 2-A) where Segev et al's theory predicts they 
should not be (point A in Fig. 1). Points C and D in Fig. 1 would correspond to coexis- 
tence of N = 2 and N = 3 domains respectively, but we find four domains and a periodic 
response (Fig. 2-C) and a quasi chaotic response (Fig. 2-D) respectively. Finally, at point 
B corresponding to the conditions proposed by SCA to illustrate their theory, one domain 
is created during each oscillation period (Fig. 2-B). In fact, the response of the system is 
very complex and within each of those zones in the Q(Iq) plane, where a constant number 
of high field domains is predicted it is possible to find all kinds of behavior: periodic with a 
single frequency, periodic with a high number of frequencies and chaotic. 

The discrepancy between the computed and predicted response should not be surprising. 
In fact, the Gunn effect is a periodic oscillation of the current (oc Jtotai) in a dc voltage 
biased semiconductor presenting negative differential velocity. It is due to periodic shedding 
and motion of charge dipole waves (high field domains) at a boundary or a nucleation site. 
It is intrinsically nonlinear so that the relevance of linear approximations such as those 
used in is questionable. We think that Kroemer's NL criterion implies that for a given 
model (equations, bias, boundary and initial conditions) it can be proved that no oscillatory 
instability is possible unless the NL product (where L is semiconductor length and N is 
doping density) is above a certain number. For Kroemer's model of the Gunn effect in 
n-GaAs under dc voltage bias see Figure 3 of M. An analytical estimate (probably not a 
very precise one) could be obtained by adapting to the Kroemer model the arguments in 
the Appendix of [|10| . This said, SCA have used a particular version of the NL criterion in 
which N is an electron density and L is either the semiconductor length (one domain) or the 
distance between domains (multiple domains). It is not surprising that this usage produces 
results that are not quantitatively correct. 

It is interesting to compare the response of the system in points C and D of Figure 
1. We observe that the dimensionless current density, j(t), in case C is periodic (Fig. 2) 
although it contains a large number of frequencies, as can be seen in the FFT spectrum 
shown in the lower right of the Figure. Instead, the response in case D is not periodic 
and the FFT spectrum appears chaotic. The difference between these two cases can be 
understood by looking at the electric field profiles as a function of time. In case C, there are 
several waves that originate at the injecting contact (right contact) and propagate towards 
the receiving contact. Each of these waves does not overtake the previous one and this 



N < N Q = 



195.565 



(15) 
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leads to a smooth, periodic response. In case D, however, a wave originates at the injecting 
contact and propaggates through the sample. Later, a new wave is injected which moves 
faster than the previous one, and that overtakes it before reaching the receiving contact. 
The interaction between waves inside the sample, leads to a complex response which may 
become chaotic. 

To understand the transition between the periodic and chaotic solutions we carried out 
numerical simulations for a fixed value of Q and varying the intensity Iq within a certain small 
range corresponding to line E in Figure 1. The results of these simulations are summarized 
in Figure 3 which shows the intensity of each mode of the FFT spectrum as a function of Iq. 
Notice that there is a frequency, which is very close to the excitation frequency, Q, which 
appears for all values of Iq. Also, all integer multiples of this frequency are present in the 
current density history. It can also be observed that there is a sharp transition between a 
periodic solution and a chaotic one. Figure 4 shows the current density history and field 
profiles for case F in the periodic region, and for case G inside the chaotic region. 

The bias condition imposes a conservation of area in the electric field profile E(z) for all 
time. Thus, the size of a travelling wave must decrease when a new wave is generated in 
the injecting contact (right contact). If the new wave grows fast enough the old wave must 
disappear to keep the area constant. In both cases, see F and G in Fig. 4, several waves are 
born at the injecting contact and propagate towards the receiving contact during a period. 
In case F, all the waves reach the left contact and this leads to a periodic response in the 
density current. As the intensity Iq increases, the wave formation velocity increases until a 
critical value above which the waves that propagate through the sample disappear before 
reaching the receiving contact. This produces a complex behavior with different patterns 
of formation-disappearance of waves. Larger values of I result in chaotic responses such as 
that shown in case G as an example. 

It is also interesting to observe that within the chaotic region there are windows in which 
the response becomes again periodic, with frecuency given by the fundamental frequency 
divided by 2, 3 or 4. In these specific ranges of Jo, the interaction between waves inside the 
sample produces a simple pattern with a periodic sequence. 

IV. CONCLUSIONS 

We have thus found that high field domains could indeed be triggered by phase-locked 
interference fringes, as suggested by SCA. However a literal use of their version of Kroemer's 
NL criterion does not often agree with numerical simulations of the model for photorefractive 
Gunn effect. In addition, we have found very interesting examples of natural and driven 
(sinusoidal interference pattern of intensity I) chaos when appropriate parameter values are 
used. In the regime for which the high field domains carry a large fraction of the bias, an 
asymptotic theory of the Gunn effect can be extended and applied to the photorefractive 
model, and used to interpret and predict the results of the numerical simulations with greater 
accuracy than SCA's linearized approach. These results will be presented elsewhere. 
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FIGURES 



FIG. 1. Level curves of intensity versus driving frequency for different integer values of Nq 
in Eq. (1). Numerical simulations have been performed for points labelled A, B, C and D, 
with dimensionless values of intensity and frecuency (io,w): A=(0. 0025, 0.05), B=(0. 035, 0.055), 
C=(0. 05, 0.0109), D=(0.07,0.05). The Fourier spectrum of the current will be displayed later in 
Figure 3 for the points on the horizontal line E. 

FIG. 2. Numerical simulations for values of the current and the oscillation frequency corre- 
sponding to points A, B, C and D in Figure 1. All other parameters correspond exactly to those 
proposed in [1] (see text). Each Figure contains four graphics in dimensionless variables: a spa- 
tio-temporal density plot of the electric field, F(y, s) (density scale on the upper right part), a plot 
of the current density j(s) (upper left), the power spectra of the j(s) signal in arbitrary units (lower 
right) and the electric field as a function of y for the fixed time marked with a white horizontal 
line in the F(y,s) density plot (lower left). 

FIG. 3. Fourier power spectra of j(t) for w = 0.05 and diferent values of iq in the range 
0.045 — 0.065, (line E of Figure 1). Each FFT is represented by a narrow horizontal band with 
gray scaled frecuency modes amplitudes: white (large) and black (small), in arbritary units. The 
arrows correspond to the simulations of a periodic signal (F) and a chaotic one (G) carried out in 
Figure 4. 

FIG. 4. Numerical simulations for the points F and G in Figure 3. The layout of the graphics 
is the same as in Figure 2. 



8 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9801256vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9801256vl 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9801256vl 



This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9801256vl 



